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Abstract 

The kinetic decoupling of dark matter (DM) from the primordial plasma sets the size of the first and smallest 
dark matter halos. Studies of the DM kinetic decoupling have hitherto mostly neglected interactions between the 
DM and the quarks in the plasma. Here we illustrate their importance using two frameworks: a version of the 
Minimal Supersymmetric Standard Model (MSSM) and an effective field theory with effective DM-quark interaction 
operators. We connect particle physics and astrophysics obtaining bounds on the smallest dark matter halo size 
from collider data and from direct dark matter search experiments. In the MSSM framework, adding DM-quark 
interactions to DM-lepton interactions more than doubles the smallest dark matter halo mass in a wide range of the 
supersymmetric parameter space. 

PACS numbers: 95.35.+d 
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I. INTRODUCTION 



The nature of dark matter is still an open question, despite the growing evidence in support of its ex- 
istence. Weakly interacting massive particles (WIMPs) are among the favorite candidates for dark matter. 
One of the unique features of WIMPs is that due to their weak interaction and heavy mass, they can lose 
thermal contact with the heat bath at a relatively early stage in the history of the Universe. The thermal, or 
as often called kinetic, decoupling of dark matter from the relativistic plasma sets the scale of the smallest 
dark matter halos, which are the first to form. Kinetic decoupling could provide a powerful cosmological 
probe on the properties of dark matter, in a way analogous to baryon decoupling, which has been unveiling 
the nature of our Universe through the baryon acoustic oscillations and the cosmic microwave background. 
Kinetic decoupling has duly received a fair amount of attention, being a compelling^ interdisciplinary av- 
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enue to explore from both particle physics and astrophysics perspectives fll- 

Kinetic decoupling is a distinct process from chemical decoupling. Kinetic equilibrium between DM 
and the primordial plasma is maintained by rapid momentum exchange through scattering. Chemical 
equilibrium, on the other hand, holds when reactions that change the number of DM particles are active 
(for WIMPs, these reactions are typically annihilation processes). Kinetic equilibrium lasts much longer 
than chemical equilibrium due to the much higher number density of relativistic plasma particles available 
for scattering compared to the non-relativistic WIMPs necessary for annihilation. At chemical decoupling 
the annihilation/creation rates fall well below the Hubble expansion rate, and the dark matter particle 
number freezes out to a constant value per comoving volume. Still the dark matter remains in thermal 
equilibrium with the relativistic plasma through frequent elastic scattering. At kinetic decoupling, even 
such scattering processes become slower than the Hubble expansion rate, and the dark matter becomes free 
from the heat bath and begins to stream freely. The dark matter kinetic decoupling determines the low- 
mass cutoff scale for the size of the dark matter halos (protohalos), which is of importance for structure 
formation in the Universe. 

In this paper, we pay particular attention to the role of DM-quark interactions in the process of kinetic 
decoupling. The role of these interactions has not been fully explored (to the extent of our knowledge, 
only Ref. Ill4ll includes them in the numerical analysis), contrary to the numerous studies that include 
scattering of DM and leptons. In addition to the lack of due attention to DM-quark interactions, another 
strong motivation for looking into the quark interactions stems from the unprecedented wealth of data 
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from the Large Hadron Collider (LHC) and direct dark matter search experiments, which provide us with 
a direct probe of the DM-quark interactions. These particle physics experiments would in principle put 
useful astrophysical bounds on the size of dark matter protohalos. 

The main results presented in this paper are (i) a kinetic-decoupling Fokker-Planck equation that ex- 
tends the expression in Bertschinger [|ll|l to general WIMP models (it was first obtained in the PhD thesis 
work of J. Kasahara under the direction of one of the authors (PG) u2]) and (ii) a connection between 
the smallest mass of the dark matter protohalos, particle searches at hadron colliders, and dark matter di- 
rect detection experiments (we illustrate this connection in the context of effective DM-quark interaction 
operators and of the Minimal Supersymmetric Standard Model). 

The layout of the paper is as follows. First, in Sec. HI we outline how we estimate the kinetic decoupling 
temperature; here we generalize the formalism presented in Ref. [11] and give a general expression for the 
momentum exchange momentum relaxation rate that can be applied to any non-relativistic WIMP. Then in 
Sec. [Ill] we apply this method to the DM-quark effective operators and show that the current collider and 
dark matter experiments set an upper bound on the smallest allowed protohalo mass. Finally in Sec. [IV] 
we quantify the importance of the DM-quark scattering relative to the DM-lepton scattering in the MSSM. 



H. THE DARK MATTER KINETIC DECOUPLING AND THE SMALLEST PROTOHALOS 

In this section we outline the formalism to estimate the kinetic decoupling temperature and the corre- 
sponding smallest dark matter protohalo mass. 

Let us first start with a heuristic order of magnitude argument for the kinetic decoupling process before 
jumping to the Boltzmann and Fokker-Planck equations. The momentum transfer per collision between 
the plasma (heat bath) at temperature T and the heavy dark matter particles % of mass m x 3> T is of order 
T, much smaller than the average momentum p of the dark matter particles, which is of order (m^T) 1 / 2 
as follows from the fact that the average kinetic energy of the dark matter particles, p 2 / (2m x ), is of order 
T. Many collisions, ~ rn^/T, are required for the dark matter to transfer a large part of its momentum 
to the plasma or to acquire it from the plasma. The momentum relaxation rate y is thus y ~ (T/m%) T e i, 
where T e i is the elastic collision rate. The dark matter is in thermal equilibrium with the plasma when the 
momentum relaxation rate is larger than the Hubble expansion rate. Thermal decoupling occurs when the 
relaxation rate becomes of the order of the Hubble expansion rate, and this defines the kinetic decoupling 
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temperature. 

We can estimate the relaxation rate and consequently the decoupling temperature more accurately 
through the Fokker-Planck equation without approximating the dark matter as a perfect fluid or fully 
collisionless gas. Ref. Ill2ll has re-derived the Fokker-Planck equation first discussed in Ref. Ill ill with- 
out assuming a specific cross section for WIMP-lepton scattering and allowing for massive particles in 
the plasma. The latter generalization is important at temperatures of the order of the electron or muon 
or quark masses. The former generalization goes beyond the zero-momentum transfer approximation of 



the previous literature [^-6 
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1411 (Mandelstam variable t = 0), and allows the study of general particle 



models. Indeed the only assumptions in Ref. [12J] are that the dark matter particle is heavy (m x 3> T 
or any other mass/energy scale) and that to momentum transfer is small. Under these assumptions, the 
Boltzmann equation for the time-dependence of the dark matter particle occupation number /x(p^) m a 
Friedmann-Robertson- Walker universe reduces to the Fokker-Planck equation 

^ -»Jr\ d (- s fi _i_ f \ , _ t 9 A 

dt 

with momentum relaxation rate 



^•S =Y(r) 4'( pz/x(i±/x)+mxr S) (1) 



Y( r) = y:^^/ ^v fi (i± fi ) P — f° dt(-t) d0x+l ^ x+l . (2) 



Here the sum extends over the species i in the relativistic plasma, with mass m ; , occupation number /, (p ; ), 
and gi statistical degrees of freedom (e.g. g q+ q = 12 for the number of spin, color, and particle-antiparticle 
states of each flavor of Standard Model (SM) quark and antiquark). The + sign in and 1 ±/j 

corresponds to bosons (stimulated emission), the - sign to fermions (Pauli blocking). And dc%+i^ x +i/dt 
is the differential scattering cross section for the elastic scattering of % and i, written as a function of the 
Mandelstam variable t and of the center-of-mass momentum p, which in the heavy % mass limit (ray 3> p) 
equals the incoming momentum of particle i in the plasma rest frame. Also, in the same limit, 

where £^ + ,^ x+ ,- is the invariant scattering amplitude and an overline indicates the usual sum over final 
polarizations and average over initial polarizations. 

The Fokker-Planck equation (OQ) automatically conserves the number of % particles per comoving vol- 
ume, since the time derivative of a 3 / x (p x ) is a total divergence in momentum space. The relaxation rate 



y(T) is an average over the thermal distribution of the mean square momentum transfer q 2 =^t times the 
collision rate. 

The expression for the momentum relaxation rate y(T) in Eq. |2]was first presented in Ref. J12I1 and it 
reproduces the formula for y(T) in Ref. Ill ill which only considered the bino scattering off massless leptons 
with a simplified invariant amplitude (the formula in Ref. IU3L 1 1411. which uses the forward scattering 
cross section do/dt\ t=0 in place of (4p 2 )~ 2 f° 4p2 dt(—t)do /dt in Eq. ©, gives a y(T) value which is 
20% larger). The formula given in Eq. |2]can be applied to a generic scattering amplitude and is thus of 
wide use. It has been implemented 112D in an extension of the DarkSUSY computer code for particle dark 
matter 1 180. 

Multiplying the Fokker- Planck equation by the % kinetic energy p^/(2m x ), and integrating in <i 3 p x 
neglecting the stimulated emission or Pauli blocking factors (1 ±/ x ~ 1), leads to an equation for the % 
kinetic temperature T%, defined as 2/3 of the average % kinetic energy, 



dt 



+ 2HT x = -2y(T) (T x -T) 



(4) 



where T is the plasma temperature. Refs. II 111 1 1311 present analytic solutions for the case y(T) proportional 
to a power of T. At temperatures greater than the kinetic decoupling temperature T^, the % particles are 
coupled to the plasma and T%~T°ca~ l . At temperatures smaller than T^, T X <*T 2 a~ 2 , as appropriate 
for non-relativistic particles of momenta p x cr x that expand freely decoupled from the rest of the uni- 
verse. In general a numerical solution of Eq. © is necessary. Here we content ourselves with estimating 



the kinetic decoupling temperature T^d as the solution of 111 , 

yiTkd) TT/rp \ 

— ~ — =n{lkd)- 



12] 



(5) 



The kinetic decoupling temperature is important to establish the mass of the smallest dark matter proto- 
halos. As pointed out in Refs. 14j,|7|,LLi|], acoustic oscillations due to the coupling between the dark matter 
and the plasma damp the amplitude of fluctuations at scales smaller than the horizon size at decoupling. 
After decoupling, dark matter particles can stream freely without interacting with the plasma, and this pro- 
cess erases fluctuations up to the distance to which they can stream from the time of kinetic decoupling. 
The mass of the smallest protohalo Mh a i ,min is determined by the larger of the DM mass inside the horizon 
at kinetic decoupling M^d and the DM mass within the free streaming length Mf s , 

Mhaio,min = max(M M , Mf s ) . (6) 
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The DM mass contained within the free streaming length Xf s is 



M f* ~ — ( T I pm0 



(7) 



where = 2%/'kf s and p m o is the dark matter density at the present time. The comoving free streaming 
scale 



rk) 

Xf s = ao dt(v/ 

Jtkd 



(8) 



where v«a 1 after kinetic decoupling, grows logarithmically during the radiation era (a t l l 2 ) and 
saturates during the matter domination era {a f 2 / 3 ). We find (see Appendix) 

TO 



where 



Okd_ = Teqf hjTeq) \ ^ 
«e 9 Tjfcrf \h{Tkd) ) 
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p eq is the total energy density at the time of matter-radiation equality. Finally, Kf s is a correction term that 
takes into account the change in the effective number of degrees of freedom between kinetic decoupling 
and the time of matter-radiation equality, 
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fTkd 



y^(Tl)h 2 / 3 (T) / , ld\nh(T) 
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where g(T) and are the energy and entropy degrees of freedom, respectively, and T* is a temperature 
much smaller than the electron-positron annihilation temperature and much larger than the temperature 
at equality, T eq < T t < 0Am e (we take T* = 1 keV for concreteness (see Appendix for the details)). For 
the dark matter velocity at decoupling Vkd, we use vm = \/6r^/5m x , with the coefficient 6/5 obtained 
numerically [11]. 



Acoustic damping is characterized by the dark matter mass inside the horizon at decoupling 
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FIG. 1: The smallest dark matter halo mass as a function the kinetic decoupling temperature Tm- The smallest 
dark halo mass is the mass within the scale characterized by the larger between the acoustic damping length and 
the free streaming length: acoustic damping at relatively small T^d, free-streaming at relatively large T^. The free 
streaming length depends on the dark matter particle mass m x through its velocity at kinetic decoupling, which scales 



as y/Tkd/m x . The feature around Tkd ~ 150 MeV is an imprint of the change in relativistic degrees of freedom during 
the QCD phase transition. 

where Tq is the present temperature, p ra o is the present DM density, and the Hubble parameter 
H(Tkd) at decoupling can be obtained from the Friedmann equation with total energy density p(7m) = 
(n 2 1 '30) g (Tkd) T* d at kinetic decoupling. 

or the effective energy and entropy degrees of freedom g(T) and h(T) we adopted the model of Ref. 



lll9ll (equation of state B) as implemented in DarkSUSY II 1811 . The factor h(T kd ) / h(To) , which was not 



present in Ref. HUH , takes into account the change of comoving volume due to the entropy increase in the 
radiation. This ratio can be bigger than 10 for temperatures of the order of the QCD scale (> 200 MeV in 
the model of our choice DIM]). 

The acoustic damping and free streaming scales are plotted in Fig. [Q For a wide range of parameters, 
the acoustic damping scale is larger than the free streaming scale and thus determines the cutoff scale of 
the smallest halo size. The free streaming length becomes more important than the acoustic damping scale 
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when Tj^/m^ becomes large as seen in the figure. This behavior is expected because the acoustic damping 
length scales as x^d ~ l/7fa*> while the free streaming length roughly scales as \fTkdfmyykd- Notice that 
the scale of the smallest dark protohalo decreases with increasing kinetic decoupling temperatures. 

HI. QUARK INTERACTIONS IN THE EFFECTIVE FIELD THEORY 

It is illustrative to consider the effective field theory approach to see the essential features of the DM- 
quark interactions and their relations to the kinetic decoupling temperature. A great feature of studying 
such effective DM-quark operators is that one can study the range of the kinetic decoupling temperature 
allowed by the recent data from the LHC and dark matter direct detection experiments that directly probe 
DM-quark interactions. 

The effective field theory approach has been studied extensively and we refer the readers to the existing 



literature for a complete survey of the effective operators [20-30] 1 . Because our purpose in this paper is to 
study the potential significance of the DM-quark scattering in the dark matter kinetic decoupling processes, 
which can be constrained from the current LHC and dark matter experiment data, we simply assume the 
dark matter is a Majorana fermion and a SM singlet, and consider the following scalar and axial-vector 
effective DM-quark point-interaction operators relevant for the direct dark matter search experiments, 

YYl 

°s = EtIxXW, (14) 

q A 

°* = E^(X/Y 5 X)(W 5 <7)- (15) 

q 

These operators lead respectively to spin-independent and spin-dependent interaction whose interaction 
strength is set by the effective cutoff scale A. The other DM-quark interaction operators besides the scalar 
and axial-vector operators vanish in the non-relativistic limit for Majorana fermion dark matter. 

We make the simplifying assumptions that the DM couples universally to the Standard Model (SM) 
up and down type quarks via 0#, and that it couples in 0$ through the quark mass suppression factor m q 



1 We do not consider the UV completion of the effective operators of our interest. Strictly speaking the effective theory breaks 
down if the 4-momentum transfer is comparable to or larger than the mass of a particle mediating the interaction (the typical 



momentum transfer is of order 7&/ for DM kinetic decoupling). See for instance Refs. [23, 30] for further discussions on the 
validity of the effective operator approach in the case of light mediators. 



8 



implied by chirality breaking. 

The spin-independent cross section per proton reads 

° SI = ^A^ (16) 
where n p = m p m%/ (m p + m x ) is the proton-DM reduced mass and / ra , e ff is the fraction of the proton mass 
m p coupled to the scalar operator O5, which is equal to the mass fraction carried by quarks plus 6/27 of 
the mass fraction carried by gluons (for the default values of the nucleon parameters in DarkSUSY, which 
we use, / OT , e ff = 0.375). The spin-dependent cross section per nucleon is 

c SD = ^nr t Aq\ 7(7+1), (17) 
where we use the value "£, q \ — 0-32 for the spin fraction carried by quarks in the nucleon [ 



32|] and 7 



is the nuclear spin (1/2 for a free nucleon). For the dark matter direct search experiment constraints, we 



use the data from SIMPLE [33] for the spin-dependent cross section and from XENON100 [13411 for the 
spin-independent interaction, which are currently the most sensitive direct detection experiments at the 
high end of the DM mass. 

A complementary set of constraints comes from collider experiments. Collider constraints do not 
suffer from the astrophysical uncertainties, such as the local galactic dark matter density and velocity 
distribution, that afflict direct search experiments. The recent LHC data release from the CMS Collabo- 



ration 113511 presented 1142 observed events in a mono jet analysis with leading jet transverse momentum 
Pt > 110 GeV, pseudorapidity |t|| <2.4 and missing transverse energy fi t > 350 GeV, to be compared 
with a Standard Model prediction of 1224 ±101 for the data sample of 4.7/fb total integrated luminosity 
at a center-of-mass energy of 7 TeV. We set a 2o collider lower bound on the effective coupling scale A. 
We implement the effective operators by treating the scalar and axial-vector operators separately (one at 
a time) in Madgraph/Madevent. We use Pythia for the hadronization and the initial/final state radiation, 
treating the jets with the pycell subroutine. 

The lower bounds on A from the direct dark matter search and collider experiments are shown in Fig. 
2 as a function of the dark matter mass m x , along with contours of the kinetic decoupling temperature T^. 
The LHC sensitivity decreases toward higher dark matter masses due to kinematic reasons. The collider 
bounds show the limits for the mono-jet events without additional nearby jets, in accord with our simpli- 
fying assumption of perfect efficiency. The direct detection experiments on the other hand suffer from the 
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(a) Scalar effective operator (b) Axial-vector effective operator 

FIG. 2: Bounds on the effective interaction scale A as a function of the dark matter mass m x for the scalar (left 
subfigure) and axial- vector (right subfigure) point-interaction effective operator. Direct detection and collider exper- 
iments exclude the regions below their respective lines (densely-dotted and dashed-dotted). These bounds become 
weaker at small and large m x respectively, due to a minimum detectable energy in direct detection experiments and 
a maximum beam energy in collider experiments. The large m x decrease of the direct detection bounds is due to 
the m~ 1 scaling of the dark matter flux onto the detector. We see that the combination of direct detection and col- 
lider bounds forces the kinetic decoupling temperature to be larger than ~ 100 MeV, a regime in which DM-quark 
scattering can be important. 



threshold to detect the DM recoil, and consequently, for the scalar operator, the collider constraints become 
stronger for light DM (m x < 10 GeV). The spin-dependent DM-nucleus scattering is not enhanced by the 
square of the atomic mass of the target nucleus, in contrast to the spin-independent interaction rate, and 
hence direct search experiments lead to bounds on the scale of the axial-vector operator that are weaker 
than the current LHC constraints . 

The combination of direct detection and collider bounds in Fig. 2 forces the kinetic decoupling tem- 
perature to be larger than ~ 100 MeV, a regime in which DM-quark scattering becomes important. The 
power-law behavior of the kinetic decoupling temperature contours in the figure emerges from the relation 
H(Tkd) = Y(7fa/)/2 and power-law dependence of the low-momentum transfer relaxation rate y(T) on A, 
T and m x . 



10 



Our effective-operator analysis provides us with upper bounds on the smallest dark protohalo mass 
directly from the current LHC and direct dark matter search limits, without having to scan the parameters 
space of a specific particle model. For instance, we find that for m x = 300 GeV the smallest allowed 
kinetic decoupling temperature is 350 MeV for the scalar operator and 150 MeV for the axial- vector 
operator, corresponding to upper limits on the smallest protohalo mass of 3 x 1O~ 9 M and 5 x 1O~ 8 M 
respectively. Note that here the protohalo cutoff mass scales as Mh a i .min 06 (Tkd^/geff(Tkd))~ 3 h e ff( Tkd ) • 
Independently of the DM mass and spin-dependent or -independent interaction, there is an absolute lower 
bounds on the kinetic decoupling temperature of ~ 100 MeV and a corresponding absolute upper bound 
of ~ 1O~ 6 M (of order of the Earth's mass) on the mass of the smallest dark protohalos. 

In the model-independent analysis of this section, each quark was decoupled at its mass scale but 
the QCD phase transition was not taken account of which is currently heavily model dependent. Indeed 
our simplified analysis adding quarks to the plasma as a free gas is not really accurate during the QCD 
phase transition, when hadrons are also present, whose treatment would require more input from lattice 
simulations. Nonetheless, our findings show that the current LHC and direct dark matter search data 
impose lower bounds on 7]y that are well in the QCD phase transition regime, supporting our statement 
that DM-quark interactions, hitherto neglected in studies of kinetic decoupling, are important and must 
be included. The inclusion of DM-quark interactions will become even more important as forthcoming 
LHC/dark matter experiments probe DM-quark interactions further, potentially pushing the lower bound 
on Tkd up and above the QCD transition regime. 



IV. QUARK INTERACTIONS IN A CONCRETE MODEL 



An advantage of specifying a concrete model rather than effective DM-quark operators is that in a 
concrete model all the DM interaction terms are present, including the interference terms, which are cum- 
bersome to specify in a model independent approach treating the different effective operators separately. 
In this section, we use a specific model to examine the relative significance of DM-quark interactions in 
comparison to DM interactions with other particles. 

As a concrete model, we choose the Minimal Supersymmetric Standard Model (MSSM), in which at 
the varying of the modelparameters the values of Tkd are known to span a wide range, from tens of MeV 
to more than a GeV [itsl [ill [l3- 17]. Because of this, the MSSM provides a good theoretical proving 
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ground to check the relative importance of the DM-quark interactions compared with the DM scattering 
off leptons that has been discussed extensively in the previous literature. 

The elastic scattering of neutralinos off the fermions in the MSSM occurs through the exchange of 
gauge bosons, Higgs bosons, and sfermions. We extended the modified DarkSUSY first developed in Ref. 
112D so as to include all possible dark matter interactions and interference terms in the MSSM for dark 
matter scattering off fermions in the primordial plasma. For the effective number of degrees of freedom, we 
adopt the model of Ref. [19] with QCD phase transition temperature of 154 MeV (equation of state B |l9|], 
the default model in DarkSUSY \iloS\ ). For the scattering of DM off quarks, the Fokker-Planck equation 
is numerically fully solved down to 154 MeV including all the relevant interactions and we simply turn 
off the DM-quark scatterings below 154 MeV. Even though the temperature at which the asymptotically 
free quark description becomes valid is expected to be higher than the QCD phase transition temperature, 
this simplification suffices for our purpose of showing the potential significance of DM-quark scattering 
in estimating the dark matter protohalo mass. Choosing a higher temperature for the threshold of DM- 
quark scattering to ensure the validity of the free quark description does not change the conclusions of 
this section (adding DM-quark scattering besides DM-lepton scattering can increase the protohalo mass 
estimation by a factor of 2 or more). 

We scanned the MSSM-7 parameter space, characterized [36] by seven parameters specified at the elec- 
troweak scale: 2 trilinear A-terms, a soft sfermion mass parameter, a gaugino mass parameter, and three 
Higgs-sector parameters. The trilinear A-terms and the soft sfermion mass matrices were assumed to be di- 
agonal to avoid flavor-changing-neutral-current issues and they were parameterized as Ay = diag(0, 0,A t ), 
Ad = rfz'ag(0,0,Ab). All the soft sfermion mass matrices are mo times the identity matrix. The gaugino 
mass GUT relation was assumed and the SU(2) gaugino mass M2 was chosen as the free parameter. The 
Higgs sector was parameterized by the CP-odd Higgs boson mass Ma, the Higgsino mass parameter /u and 
the ratio tan [3 of Higgs vacuum expectation values. These seven parameters were randomly scanned over 
the ranges 10 < tan [3 < 50, < \A t j>\ < 10 TeV and [lOGeV, 107W] for the remaining mass parameters 
mo, Mi, Ma, and /j. We applied the phenomenological bounds available in the DarkSUSY subroutines II 1 8f1 
for the theoretical consistency of the model and the experimental constraints. We excluded parameter sets 
that have charged or colored vacua, that have potentials unbounded from below, that do not have the neu- 
tralino as the lightest supersymmetric particle, and that violate experimental bounds on supersymmetric 
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masses or rare processes such as b — > sy. 

Because our purpose here is to show quantitatively the significance of the DM-quark scattering we 
do not require that the neutralino thermal relic abundance should match the observed cold dark matter 
value (actually, through non-thermal production, entropy production, or non-standard expansion history, 
the neutralino thermal relic abundance can be smaller or larger than the cosmological value Q-dm ~ 0.1, 
even if the neutralinos comprise all of the cold dark matter Q37|]). 
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FIG. 3: The neutralino kinetic decoupling temperature as a function of the neutralino mass for a random scan in the 
MSSM-7 parameter space. The inclusion of DM-quark scattering can change the kinetic decoupling temperature by 
30% and hence the dark matter protohalo mass by more than a factor of 2. 

In Fig. [3] we compare the neutralino kinetic decoupling temperatures with and without including DM- 
quark interactions for a random scan in the MSSM-7 parameter space. The open circles represent the 
Tkd values computed including the DM-quark interactions, and the crosses are the Tm values computed 
without including the DM-quark interactions. We first observe that kinetic decoupling can occur above the 
QCD phase transition scale, as first pointed out in Ref. (SO, which performed an analogous numerical scan 
without including the DM-quark interactions. For < Tqcd, the quarks are bounded inside the hadrons 
and the pions were much less abundant than the light leptons in the thermal bath. We numerically checked 
that the inclusion of pion scattering affects the kinetic decoupling temperature by less than one per cent 
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in our scenarios 2 . Fig. [3] shows that can be significantly influenced by DM-quark interactions for 
Tkd above the QCD scale. In our MSSM-7 parameter scan we found that the typical ratio of T^d with and 
without DM-quark interactions ranged between 1 and ~ 1.3. When the smallest protohalo mass Mhai ,min 
is set by the acoustic damping scale, it scales with as Mh a io,min ~ T M . Hence the inclusion of DM- 
quark scattering can lead to a factor 2 or more increase in the smallest protohalo mass for a wide range of 
MSSM parameter values. When the smallest protohalo mass Mh a io,min is set by the free streaming length, 
it roughly scales as Mh a io,min ~ T kd ' , and it less affected by the inclusion of DM-quark interactions. 

This conclusion on the importance of including DM-quark interactions would not change even if we 
had switched off DM-quark scattering up to temperatures four or five times higher than the QCD phase 
transition temperature, as to guarantee the presence of free quarks, because our MSSM-7 parameter scan 
already covers such regime of 7^ well above the QCD phase transition temperature 3 . 



V. DISCUSSION/CONCLUSION 



Studies of particle dark matter can provide us with a unique link between particle physics and astro- 
physics. To illustrate one such connection, we examined the importance of DM-quark interactions in the 
kinetic decoupling of particle dark matter and the consequent mass of the smallest dark protohalos. We 
focus on a model independent analysis through effective DM-quark interaction operators and a model 
specific analysis within a seven-parameter MSSM. 

In the effective operator approach, in which we assumed DM-lepton interactions are negligible, we 
found that current direct dark matter search and collider constraints force the DM kinetic decoupling 
temperature to exceed ~ 100 MeV, whether the interaction is spin-dependent or spin- independent. 



1 The leading order DM-pion coupling term cn XX can be estimated with a coupling constant [9 

i 

Tfl 

c= T( * A fu+fd), (18) 
2{m u + m d ) 

where /„, with q — up or down quark in this case, is the coefficient of the effective scalar current operator f q qq%X in the 
MSSM 1 10]. This form of cn 2 %% was derived using the soft-pion technique in the framework of chiral perturbation theory, 
which is valid for an energy scale < 4%/% (f K is the pion decay constant), and hence is applicable to the soft-pion case. 
3 That such high Tkd are common can be seen for instance in the MSSM scans in Ref. 11411 . which used a different QCD phase 
transition model and a different method of estimating the kinetic decoupling temperature. 
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This sets an upper limit of ~ 1O -6 M (of order of the Earth's mass) to the mass of the smallest dark 
protohalo. 

In the MSSM-7 study, we found that inclusion of neutralino-quark interactions can increase the smallest 
dark protohalo mass by more than a factor of 2 whenever the neutralino decoupling temperature exceeds 
the QCD transition temperature and acoustic damping dominates over free streaming. 

The constraints on the kinetic decoupling temperature from the current LHC and direct dark matter 
search experiments turned out to be in the regime of the QCD phase transition. The forthcoming LHC and 
direct dark matter search data will most likely push this 7Jy bound up, possibly beyond the quark-hadron 
transition, into a regime in which our approximation of a gas of free quarks is definitely applicable. 

If DM-quark interactions are discovered at the LHC or in direct detection experiments, measurements 
of their strength can provide us with a lower bound of the mass of the smallest dark matter protohalos, 
once proper account is taken of possible additional interactions with leptons. 

The survival of the smallest dark matter protohalos to the present time, and their observability, depend 
strongly on complicated tidal forces and other astrophysics such as stellar interactions, and have been the 
subject of vigorous debate. More input from computer simulations of structure formation in the universe 
would help clarify these issues, despite the vast dynamic range of masses and non-linear effects involved 
in simulating the relevant processes 

Although we discussed only the hadron collider (LHC) data that constrain the DM-quark interaction 
operators, an analogous exercise can be performed for the DM-lepton interactions considering a lepton 



collider Q38 



390 which is crucial for a concrete model such as the MSSM where DM-lepton interactions 



cannot be neglected for the estimation of the kinetic decoupling. Inclusion of non-collider experiments 
and indirect search experiments also deserve further studies. For instance the constraints from the cosmic 

nn 

antiproton flux can put relatively tight constraints on the axial- vector operator |J40L 14 1 1] - and more system- 
atic studies of the dark matter kinetic decoupling including additional experiments and additional effective 
operators beyond those considered here are worthy of future work. 

Our results on the impact of collider and dark matter search experiments onto the formation of structure 
in the early universe calls for further exploration of the creation and evolution of dark matter protohalos, 
in view of their potential role in probing the nature of the dark matter. 
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Appendix A: Derivation of the free streaming length formula 

Here we sketch the derivation of Eq. © for the free streaming length, which includes the change in 
the effective degrees of freedom before matter-radiation equality. Since the expressions for the energy and 
entropy degrees of freedom g(T) and h(T) usually do not extend down to temperatures ~ T eq , we divide 
the integration in 

ft* v 

Xf s = a dt- (Al) 

Jt kd a 

into two parts, which we join at a temperature T* defined so that aT = const for T <T*. The tempera- 
ture r* is thus much smaller than the electron-positron annihilation temperature but much larger than the 
temperature at equality. 

For T < I 7 *, we use an FRW model with matter and radiation, 




(A2) 



where H eq = ^SnGp eq /3 is the Hubble parameter at the matter-radiation equality (a = a eq ), with p eq 
equal to the total density at that time (the contribution to X/s from the low-redshift cosmological constant 
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term is negligible). We find 



r° v no dt 

/ dt - = aov kd a kd / -j 
Jn a J u a L 

f a o da 
= aov kd a kd / ~ 

Ja t a 3 H(a) 

- y/2 a ° Vkdakd ' 



a liqHeq 



__ VkdOkd 

where, on using aT = const for a* < a < ao, 

Cx = 

with 



In- 



In 



Cx 



2 + Cx 



yT+q-1 
VT+a+l 

to 

> 



a=a /a e , 



a=a*/a e , 



(A3) 
(A4) 

(A5) 

(A6) 



7 

1 + ^-1. 

r 



(A7) 



, . (A8) 
aovo flo v j 

For r > r*, we include the full temperature dependence of the degrees of freedom. Using conservation 
of total entropy a 3 T 3 h(T) = const, from which 



da 



( _ldlnh(T)\ dT 



3 dlnT 



(A9) 



we find 



a 



("d, v - = 

Jt kd a 



ttQV kd a kd 


ra* 


da 


'akd 


a i H(a) 


aov kd a kd 




a 2 H{a*) da 


alH{a*) 


Ja kd 


a 2 H(a) a 


v kd a kd 


1 


y/gfr) h 2 '\T) 


2a eq C / 


1 + 


2p h h 2 l 3 {T,) y/gp) 



(A10) 
(All) 
(A12) 



In the last step, we wrote H(T*) in terms of H eq , and then C, using the matter and radiation model at T < T* 
where aT = const, 



(A13) 
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After adding Eqs. (IA61 ) and (|A121 ). we formally extend the integration in Eq. (IA61 ) from x* to using 
the definition 



CX: 



T 

l + i£i_l 

r 



(A14) 



and the mathematical identity 







fin CX " 




2 + Cx_ 





l kd 



1 dT 



T T 



(A15) 



Notice that since dT ^ const for T > the variable x in Eq. (IA14I) is not the conformal time at plasma 
temperatures T > T*, as instead it is at T < T*, but is only a convenient mathematical variable. This gives 



■fs 



In- 



Cx 



2 + Cx 



to 



(A16) 



with 



K 



fs 



1 + 



= f T " Vg(Tdh 2 /*(T) f U 



d]nh(T)\ dT f T kd 1 



dlnT 



dT 



T T 
1 + ^ 



(A17) 



The approximation ^> T eq then gives the expression of Kt s in Eq. (fTTI) . 
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